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Abstract — The  main  purpose  of  this  work  is  to  investigate 
the  porous  medium  and  diffusion  thermo  effects  on 
unsteady  combined  convection  magneto  hydrodynamics 
boundary  layer  flow  of  viscous  electrically  conducting 
fluid  over  a vertical  porous  surface  in  the  presence  of  first 
order  chemical  reaction  and  thermal  radiation.  The  slip 
boundary  condition  is  applied  at  the  porous  surface.  A 
uniform  Magnetic  field  is  applied  normal  to  the  direction 
of  the  fluid  flow.  The  non-linear  coupled  partial 
differential  equation  are  solved  by  perturbation  method 
and  obtained  the  expressions  for  concentration, 
temperature  and  velocity  fields.  The  rate  of  mass  transfer 
in  terms  of  Sherwood  number  h , the  rate  of  heat  transfer 
in  terms  of  Nusselt  number  Nu  and  the  Skin  friction 
coefficient  Cf  are  also  derived.  The  Profiles  of  fluid  flow 
and  derived  quantities  for  various  values  of  physical 
parameters  are  presented  and  analyzed. 

Keywords — Diffusion-thermo  effect,  Thermal  radiation, 
Chemical  reaction,  Magnetic  field,  Absorption  of 
Radiation. 

Nomenclature 

A Suction  velocity  parameter 

B0  Magnetic  field  of  uniform  strength 

C Species  concentration 

C*  Dimensional  concentration 

Cp  Specific  heat  at  constant  pressure 

Cs  Concentration  susceptibility 

Cw  Species  concentration  at  plate 

^ Species  concentration  far  away  from  the  plate 
Cf  Skin  friction  coefficient 

D Molecular  diffusivity 

ebX  Planck’s  function 

g Acceleration  due  to  gravity 


Gm  Solutal  Grashof  number 

Gr  Thermal  Grashoff  number 

Dm  Coefficient  of  mass  diffusivity 

Kt  Thermal  diffusion  ratio 

Kxw  Absorption  coefficient  at  the  wall 

M Magnetic  field  parameter 

F Radiation  parameter 

Nu  Nusselt  number 

n*  Constant 

P*  Pressure 

Pr  Prandtl  number 

Q Heat  source/sink  parameter 

Q*  Sink  strength 

Q0  Dimensional  heat  absorption  coefficient 

Q1  Absorption  of  radiation  parameter 

Q j Coefficient  of  proportionality  for  the  absorption 

q Radiative  heat  flux 

Sc  Schmidt  number 

Sh  Sherwood  number 

T Fluid  Temperature 

T*  Temperature  of  the  fluid  near  the  plate 

Tw  Temperature  at  the  wall 

Tm  Temperature  far  away  from  the  plate 

u*  , v*  Components  of  dimensional  velocities 

u*  Velocity  of  the  fluid  along  x* 

Uf  Free  stream  dimensional  velocity 

Up  W all  dimensional  velocity 

v*  Velocity  of  the  fluid  along  y* 

U0  Scale  of  free  stream  velocity 

V0  Section  velocity 

Cp  Specific  heat  at  a constant  pressure 

Cv  Specific  heat  at  a constant  volume 

x* , y*  Dimensional  the  distances  along  and 

perpendicular  to  the  plate 
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Rex  Local  Reynolds  number 

Greek  Ssymbols 

p Fluid  dynamic  viscosity 

a Fluid  thermal  diffusivity 

v Coefficient  of  kinematic  viscosity 

P Coefficient  of  Volumetric  expansion  for  the  heat 

transfer 

/?*  Coefficient  of  Volumetric  expansion  for  the  fluid 
/?„  Magnetic  field  coefficient 

(ir  , Pc  Thermal  and  concentration  expansion  coefficients 
p Density  of  the  fluid 

a Electrically  conductivity  of  the  fluid 

rj  Dimensionless  normal  distance 

Fluid  currents  formed  in  a fluid-saturated  porous  medium 
during  convective  heat  transfer  have  many  important 
applications,  such  as  oil  and  gas  production,  central  grain 
storage,  porous  insulation,  and  geothermal  energy.  The 
study  of  natural  convection  through  porous  medium  also 
throws  some  light  on  the  influence  of  environment  such  as 
temperature  and  pressure  on  the  germination  of  seeds.  In 
many  situations,  the  heat  and  mass  transfer  on  the 
hydromagnetic  flow  near  the  vertical  plate  is  encountered 
e.g.,  the  cooling  of  nuclear  reactor  with  electrically 
conducting  coolants  such  as  liquid  sodium  and  mercury, 
studied  by  Rath  and  Parida  [1].  Raptis  [2]  , Jha  and  Prasad 
[3]  have  studied  the  steady  free-convection  flow  and  mass 
transfer  through  a porous  mediumbounded  by  an  infinite 
vertical  plate  for  the  flow  near  plateby  using  the  model  of 
Yamamoto  and  Iwamura  [4],  Gebhart  and  Pera[5]  have 
studied  the  laminar  flows  which  arise  in  fluids  due  to  the 
interaction  of  the  force  of  gravity  and  density  differences 
caused  by  the  simultaneous  diffusion  of  thermal  energy 
and  of  chemical  species. 

As  stated  by  Pal  and  Talukdar  [6],  convection  in  porous 
media  has  gained  significant  attention  in  recent  years 
because  of  its  importance  in  engineering  applications  such 
as  solid  matrix  heat  exchangers,  geothermal  systems, 
thermal  insulations,  oil  extraction  and  store  of  nuclear 
waste  materials.  Convection  in  porous  media  also  is 
applied  to  underground  coal  gasification,  ground  water 
hydrology,  iron  blast  furnaces,  wall  cooled  catalytic 
reactors,  solar  power  collectors,  energy  efficient  drying 
processes,  cooling  of  electronic  equipment’s  and  natural 
convection  in  earth’s  crust.  Reviews  of  the  applications 
associated  to  convective  flows  in  porous  media  can  be 
found  in  Nield  and  Bejan  [7]. The  fundamental  problems 
of  flow  through  and  past  porous  media  has  been  studied 
extensively  over  the  years  both  theoretically  and 
experimentally  by  Cheng  [8]  and  Rudraiah[9].The  effect 
of  Beavers-Joseph  slip  velocity  and  transverse  magnetic 
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0 Heat  source  parameter 

0X  Slip  parameter 

k Thermal  conductivity 

Supersripts 

Differentiation  with  respect  to  y 
* Dimensional  properties 

Subscripts 

p Plate 

w W all  condition 

oo  Free  stream  condition 

I.  INTRODUCTION 

field  on  an  electrically  conducting  viscous  fluid  in  a 
horizontal  channel  bounded  on  both  sides  porous 
substrates  of  finite  thickness  was  studied  Rudraiah  et  al. 
[10],  which  is  equivalent  to  the  problem  of  forced 

convection  where  the  momentum  equation  is  independent 
of  concentration  distribution  and  the  diffusion  equation  is 
coupled  with  the  velocity  distribution  using  Beavers- 
Joseph  slip  condition  at  the  porous  interface. 

In  many  chemical  engineering  processes,  the  chemical 
reaction  occurs  between  a foreign  mass  and  the  fluid  in 
which  the  plate  is  moving.  These  processes  take  place  in 
numerous  industrial  applications  e.g.,  manufacturing  of 
ceramics  or  glassware,  polymer  production,  and  food 
processing.  Chemical  reactions  can  be  codified  as  either 
homogeneous  or  heterogeneous  processes.  This  depends 
on  whether  these  occur  at  an  interface  or  as  a single  phase 
volume  reaction.  Many  transport  processes  exist  in  a 
nature  and  mass  transfer  as  a result  of  combined  buoyancy 
effects  of  thermal  diffusion  and  diffusion  of  chemical 
species.  Apelblat  [11]  studied  analytical  solution  for  mass 
transfer  with  a chemical  reaction  of  the  first  order.  Das  et 
al.  [12]  have  analyzed  the  effects  of  homogeneous  first 
order  chemical  reaction  on  the  flow  past  an  impulsively 
started  infinite  vertical  plate  with  uniform  heat  flux  and 
mass  transfer.  Chambre  and  Young  [13]  have  analyzed  a 
first  order  chemical  reaction  in  the  neighborhood  of  a 
horizontal  plate. 

Radiative  convective  flows  are  encountered  in  various 
ways  in  the  environment  e.g.,  heating  and  cooling 
chambers,  evaporation  from  large  open  water  reservoirs, 
fossil  fuel  combustion  energy  processes,  solar  power 
technology,  astrophysical  flows,  and  space  vehicle  re- 
entry. Radiative  heat  and  mass  transfer  play  an  important 
role  in  manufacturing  industries  for  the  design  of  reliable 
equipment.  Nuclear  power  plants,  gas  turbines,  and 
various  propulsion  devices  for  satellites,  missiles,  aircraft 
and  various  space  vehicles  are  examples  of  such 

Page  I 66 


International  Journal  of  Advanced  Engineering  Research  and 
https://dx.doi.ors/10.22161/iiaers. 4. 5. 12 
engineering  applications.  If  the  temperature  of  the 
surrounding  fluid  is  rather  high,  radiation  effects  play  an 
important  role  and  this  situation  does  exist  in  space 
technology.  In  such  cases,  one  has  to  take  into  account  the 
effect  of  thermal  radiation  and  mass  diffusion. 

The  thermal  radiation  effects  of  an  optically  thin  gray  gas 
bounded  by  a stationary  vertical  plate  was  studied  by 
England  and  Emery  [14].  The  radiative  natural  convective 
flow  of  an  optically  thin  gray-gas  past  a semi-infinite 
vertical  plate  was  considered  by  Soundalgekar  and  Takhar 
[15],  In  all  above  studies,  the  stationary  vertical  plate  is 
considered.  The  effects  of  the  thermal  radiation  and  free 
convection  flow  past  a moving  vertical  plate  studied  by 
Raptis  and  Perdikis  [16].Ibrahim  et  al.  [18]have  recently 
reported  computational  solutions  for  transient  reactive 
magneto  hydrodynamic  heat  transfer  with  heat  source  and 
wall  flux  effects.  They  have  also  analyzed  the  effects  of 
the  chemical  reaction  and  radiation  absorption  on  the 
unsteady  MHD  free  convection  flow  past  a semi-infinite 
vertical  permeable  moving  plate  with  heat  source  and 
suction.  Due  to  the  importance  of  Soret  (thermal- 
diffusion)  and  Dufour  (diffusion  thermo)  effects  for  the 
fluids  with  very  light  molecular  weight  as  well  as  medium 
molecular  weight  many  investigators  have  studied  and 
reported  results  for  these  flows  of  whom  the  names  are 
Dursunkaya  and  Worek  [19],  Anghel  et  al.  [20], 
Postelnicu  [21]  are  worth  mentioning. 

Chamkha[24]  investigated  unsteady  convective  heat  and 
mass  transfer  past  a semi-infinite  porous  moving  plate 
with  heat  absorption.  Chamkha[25]  studied  the  MHD  flow 
of  a numerical  of  uniformly  stretched  vertical  permeable 
surface  in  the  presence  of  heat  generation/absorption  and  a 
chemical  reaction.  Mohamed  [26]  has  discussed  double 
diffusive  convection  radiation  interaction  on  unsteady 
MHD  flow  over  a vertical  moving  porous  plate  with  heat 
generation  and  Soret  effects.  Muthucumaraswamy  and 
Janakiraman[27]  studied  MHD  and  radiation  effects  on 
moving  isothermal  vertical  plate  with  variable  mass 
diffusion. Rajesh  and  Varma[28]  studied  thermal  diffusion 
and  radiation  effects  on  MHD  flow  past  a vertical  plate 
with  variable  temperature. Kumar  and  Varma[29] 
investigated  thermal  radiation  and  mass  transfer  effects  on 
MHD  flow  past  an  impulsively  started  exponentially 
accelerated  vertical  plate  with  variable  temperature  and 
mass  diffusion.  Raptis  et  al.[30]  studied  the 
hydromagnetic  free  convection  flow  of  an  optically  thin 
gray  gas  taking  into  account  the  induced  magenetic  field 
in  the  presence  of  radiation  and  the  analytical  solutions 
were  obtained  by  perturbation  technique.  Orhan  and 
Ahmad  [31]  examined  MHD  mixed  convective  heat 
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transfer  along  a permeable  vertical  infinite  plate  in  the 
presence  of  radiation  and  solutions  are  derived  using 
Kellar  box  scheme  and  accurate  finite -difference  scheme. 
Ahmed  [32]  investigated  the  study  of  influence  of  thermal 
radiation  and  magnetic  Prandtl  number  on  the  steady 
MHD  heat  and  mass  transfer  mixed  convection  flow  of  a 
viscous,  incompressible,  electrically-conducting 
Newtonian  fluid  over  a vertical  porous  plate  with  induced 
magnetic  field. 

When  heat  and  mass  transfer  occur  simultaneously  in 
moving  fluid,  the  relations  between  the  fluxes  and  the 
driving  potentials  are  of  a more  intricate  nature.  It  has 
been  observed  that  an  energy  flux  can  be  generated  not 
only  by  temperature  gradients  but  also  by  concentration 
gradients.  The  energy  flux  caused  by  a concentration 
gradient  is  termed  the  diffusion-thermo  (Dufour)  effect. In 
most  of  the  studies  related  to  heat  and  mass  transfer 
process,  Soret  and  Dufour  effects  are  neglected  on  the 
basis  that  they  are  of  a smaller  order  of  magnitude  than  the 
effects  described  by  Fourier’s  and  Fick’s  laws.  But  these 
effects  are  considered  as  second  order  phenomena  and 
may  become  significant  in  areas  such  as  hydrology, 
petrology,  geosciences,  etc.  For  fluids  with  medium 
molecular  weight  (H2,  air),  Dufour  and  Soret  effects 
should  not  be  neglected  as  indicated  by  Eckert  and  Drake 
[17]  Prakash  etal.  [33]  analyzed  Diffusion-Thermo  and 
Radiation  Effects  on  Unsteady  MHD  Flow  through  Porous 
Medium  Past  an  Impulsively  Started  Infinite  Vertical  Plate 
with  Variable  Temperature  and  Mass  Diffusion.  To  our 
best  knowledge,  the  interaction  between  the  diffusion- 
thermo  and  chemical  reaction  in  the  presence  of  porous 
medium,  heat  absorption/generation,  thermal  radiation, 
radiation  absorption  effects  has  received  little  attention. 
Hence,  an  attempt  is  made  to  study  the  diffusion-thermo 
effects  on  an  unsteady  MHD  free  convective  heat  and 
mass  transfer  flow  of  a viscous  incompressible  electrically 
conducting  fluid  through  porous  medium  from  a vertical 
porous  plate  with  varying  suction  velocity  in  slip  flow 
regime. 

To  the  best  of  our  knowledge,  the  interaction  between  the 
diffusion-thermo  and  chemical  reaction  in  the  presence  of 
porous  medium,  heat  absorption/generation,  thermal 
radiation,  radiation  absorption  effects  has  received  little 
attention.  Hence,  an  attempt  is  made  to  study  the 
diffusion-thermo  effects  on  an  unsteady  MHD  free 
convective  heat  and  mass  transfer  flow  of  a viscous 
incompressible  electrically  conducting  fluid  through 
porous  medium  from  a vertical  porous  plate  with  varying 
suction  velocity  in  slip  flow  regime. 
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II.  MATHEMATICAL  FORMULATION 

Consider  an  unsteady  two  dimensional  flow  of  an 
incompressible  viscous,  electrically  conducting  and  heat- 
absorbing fluid  past  a semi-infinite  vertical  permeable 
plate  embedded  in  a uniform  porous  medium  and 
subjected  to  a uniform  transverse  magnetic  field  in  the 
presence  of  thermal  and  concentration  buoyancy  effects. 
The  applied  magnetic  field  is  also  taken  as  being  weak  so 
that  Hall  and  ion  slip  effects  may  be  neglected.  We 
assume  that  the  Dufour  effects  may  be  described  by  a 
second-order  concentration  derivative  with  respect  to  the 
transverse  coordinate  in  the  energy  equation.  Further  to 

(1) 


our  assumption  that  there  is  no  applied  voltage  which 
implies  the  absence  of  an  electric  field.  The  plate  is 
maintained  at  constant  temperature  Tw  and 
concentrationCw,  higher  than  the  ambient  temperature  Tm 
and  concentration  C <*,  respectively.  The  chemical  reactions 
are  taking  place  in  the  flow  and  all  thermo  physical 
properties  areassumed  to  be  constant.  Due  to  the  semi- 
infinite plane  surface  assumption,  the  flow 
variables  are  functions  of  y*and  time  t*  only. Under  the 
usual  boundary  layer  approximations  the  governing 
equations  are  governed  by  the  following  equations. 


(2) 


f,  = 0 

dy* 

du*  * du*  1 dp*  d2u*  oBq  * n f . n f . 

JF  + V -^  = --p^  + v^-—u  +^r(.T  — Too)  + gPc(c  -c„) 

S + = — S'— (7’*-7’oo)  + ^(C*-Coo)+^^^(3) 

at-  dy * pCpdy'2  pCpdy * pCpy  007  PCp  v 007  Cs  pCp  Sy*2v 

dC*  dC*  d2  C* 

— + v*  — — D — — pi  — R(C*  — Cqo)  (4)wherex*,  y*  and  t*are  the  dimensional  distances  along  x*  and  y*directions 

and  dimensional  time  respectively,  u*  and  v*  are  the  components  of  dimensional  velocities  along  x*  and 
y*directionsrespectively,  T*is  the  dimensional  temperature,  C*is  the  dimensional  concentration,  Tw  and  Cware  the  temperature 
and  concentration  at  the  wall,  C^and  Tm are  the  free  stream  dimensional  concentration  and  temperature,  p is  the  density,  v is 
kinematic  viscosity,  Cp  is  the  specific  heat  at  constant  pressure,  a is  the  fluid  electrical  conductivity,  B0  is  the  magnetic 
induction,  If  * is  the  permeability  of  the  porous  medium,  q)  is  radiative  heat  flux,  Q0  is  the  dimensional  heat  absorption 
coefficient,  Qiis  the  coefficient  of  proportionality  for  the  absorption,  R is  the  chemical  reaction,  j3T and  /?care  the  thermal  and 
concentration  expansion  coefficients,  D is  the  molecular  diffusivity,  Dm is  the  coefficient  of  mass  diffusivity,  KT is  the  thermal 
diffusion  ratio, 

Q0(T*  — 7’00)is  assumed  to  be  the  amount  of  heat  generated  or  absorbed  per  unit  volume  and  Q0  is  a constant. 

The  radiative  heat  flux  is  considered,  which  is  given  by  Cogley  et  al.  [22],  Pal  and  Talukdar  [6]  as 


= 4(7*  - Tay 


(5) 


where/'  = f()  KXw  dA,  /f;iVVthe  coefficient  of  absorption  near  the  wall  and  ehx  is  Planck’s  function. 

Under  the  above  stated  assumption,  the  initial  and  boundary  conditions  for  the  velocity  distribution  involving  slip  flow, 
temperature  and  concentration  distributions  are  defined  as: 

VfcdU*  rr,  , _rrr,  rf,  N n«t« 


u — u 


slip  ~ - Tw  + e(Tw  ~ T^e1 


aty*  = 0 (6) 

C*  = Cw  + e(Cw  - C00)en‘t\ at y*  = 0 (7) 

u*  = //^  = //0(l  + een*t‘),r*->7OTi 
C*->C00,asy‘->ra  (8) 

From  Eq.  (1)  it  is  clear  that  the  suction  velocity  at  the  plate  surface  is  either  constant  or  a function  of  time  only.  Hence,  it  is 
assumed  that 

d*  — — Vq(1  + eAe11*1*)  (9) 

where  F(lis  the  mean  suction  velocity  and  eA  « 1.  The  negative  sign  indicated  that  the  suction  velocity  is  directed  towards  the 
plate. 

In  the  free  stream  Eq.(2)  gives 

1 dp* 


“-F  tA-Upp  — ^2  JJ* 

p dx * dt*  p 0 00 


(10) 


On  introducing  the  non-dimensional  quantities 


u*  D*  v 0y* 

u = — v = — ,y  = 

U0  v0’y  v 


u =U-± 

’ 00  a ’ 
uo 
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_ vft*  r - To,,  _ c*  — ^ 

t — , u — ,L  — 


T —T 

1 w 1 cc 


r — r * 

'-'W  00 


n*v  pgv(Tw-Tm)pr 

n = — , Gr  = ^ , 

V02  « 

pgv(Cw  - C^Pc  avBg 
Gm  ,M  — 


vCv 

Pr  = ~~ , 0 = 

K 


GqVq 

QqV 

pv02cp 


,Qi 


F = 


4 vl'  v 

— ~ ,Sc  — — ,v 

VgPGp  D’y 


pv02’ 

QiV(Gw  - Cqq) 
V02(Tw~Tm)  ' 
Rv 


V02’ 


Du  = 


Dm^-T  (Cw  Cqo) 


(ii) 


CSK  (Tw  Too) 

In  view  of  the  above  dimensionless  variables,  the  basic  field  of  Eqs.  (2)  to  (4)  can  be  expressed  in  dimensionless  form  as 


^ - (1  + eAent) P = — + M(f/oc  -U)  + ^-+Gr9  + GmC, 


at 

ae 

at 


dy 

dy 


dt 
1 d20 
Pr  dy2 


dy2 


^-(l  + eAent)^  = ^2-F6  + Q1C -09  + ^0 


(12) 

(13) 


^-(l  + eAent)^  = ^-yC(14) 


dt  dy  Sc  dy2 

The  corresponding  initial  and  boundary  conditions  in  equations  (6)-(8)  in  non-dimensional  form  are  given  below 

nnt 


U — U 


du 


slip  — gy ' 9 — 1 A ee 


C — 1 + eentaty  = 0 
asy  -*  oo 


(15) 

U Coo  = (1  + eent)  ,9->0,C 

(16) 


III.  METHOD  OF  SOLUTION 

The  equations  (12)  to  (14)  are  coupled  non-linear  partial  differential  equations  whose  solutions  in  closed-form  are  difficult  to 
obtain. To  solve  these  coupled  non-linear  partial  differential  equations,  we  assume  that  the  unsteady  flow  is  superimposed  on 
the  mean  steady  flow,  so  that  in  the  neighbourhood  of  the  plate,  we  have 

u = /o(y)+een7i(y)  + 0(e2)  (17) 

o = 5o(y)+e  entgt(y)  + 0(£2)  (18) 

C = h0(y)+E  entht(y)  + 0(E2)  (19) 

By  substituting  the  set  of  Eqs.  ( 17)-(  19)  into  Eqs.(12)-(14)  and  equating  the  harmonic  and  non-harmonic  terms,  and  neglecting 

the  higher  order  terms  in  £,  we  obtain 

/o"  + /o'  - MJ0  = -Mx  - Grg0  - Gmh0  (20) 

fi  + A - M/i  = -(M  + n) 

—Afg  — Grg1  — Gmh±  (21) 

g'o  + Prg'o  - Pr (F  + 0)go  = -PrQ-^hg  - Duh'd  (22) 

g'l  + Prg[  — Pr(F  + 0 4-  ri)gx  = —PrQ1h1  — Duh'd  — PrAg'0 
(23) 

h'd  + S ch'g  - Scyhg  = 0 (24) 

ft"  + Sch[  — Sc(y  + n )h1  — —ASch'0(25) 

where  the  prime  denotes  the  differentiation  with  respect  to  y.  Now  the  corresponding  boundary  conditions  are 

fo  — 0i/o</i  = 0i/i<5o  = 1> 

g1  = 1,  h0  = 1,  = 1 aty  — 0 (26) 

fo  = l.A  = l-5o  0,9i  -»  0,ho  ->  0, 

->  0,  as  y ->  oo  . (27) 

By  solving  the  set  of  Eqs. (20)-  (25)  with  the  help  of  boundary  conditions  (26)-(27),  the  following  set  of  solutions  are 
obtained 

fo  = 1 + B4e~m<>y  - A3e-m?y  - 42e~p^( 28) 
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ft  = 1 + B5e~miy  + A4e~m6y  - Awe~m 2y  - B1e~Piy  - B2e~m ^ - B3e~m *y  (29) 
g0  = (l-Af)e-m2y + Aie-p^m 

g4  = Age-msy  + A6e~piy  + A7e~m2y  + Aae~m*y  (31) 

h0  = e~Piy  (32) 

hx  = (1  - A5)e~m*y  + A5e~Piy  (33) 
u(y,  t)  = 1 + B4e-m*y  - A3e~m2y  - A2e~p^y+E  ent(l  + B5e~m2y  + A 4e-m^y  - Awe-m2y  - Bxe~p'y  - B2e~m^y  - 
B3e~m*y)  (34) 

B(y,  t)  = (1  - A1)e~rri2y  + Aie-piy+ 

£ ent(Age~m5y  + A6e~piy  + A7e~m2y  + 48e“m^)  (35) 

C(y,  t)  = e-Piy+£  ent((l  — A5)e~miy  + i45e~Piy)  (36) 


The  coefficient  of  Skin-friction,  the  rate  of  heat  transfer  in  terms  of  Nusselt  number  and  the  rate  of  mass  transfer  are  the 
important  physical  parameters  for  this  kind  of  boundary  layer  flow.  Hence  these  quantities  are  defined  and  derived  as  follows : 


Cfx  — 


pu0v0 


= (t;)  = ~B4m6  + A3™2  + A2Pi+£  ent(— f?5m3  - A4m6  + Awm2  + B1P1  + B2m5  + B3m4 ) 

Koy/ aty= 0 


(37) 


Nur  = x 


Nux j =(^\  =(-m2(l-A1)-A1P1)+Eent(-m5Ag-P1A6-m2A7-m4Aa) 

f\Cx  \Oy/  nt-w-n 


(Tw  Too) 


aty= 0 


(38) 


= ~Pi+  G ent(—m4(l  - As)  - AsPf)  (39) 
Where!?  ex  = — is  the  local  Reynolds  number. 


Shx 


ydy’Jaty= 0 

* (Cw  - COT)’ 


aty= 0 


IY.  RESULTS  AND  DISCUSSION 

The  analytical  solutions  are  performed  for  concentration, 
temperature  and  velocity  for  various  values  of  fluid  flow 
parameters  such  as  Schmidt  number  Sc  , chemical 
reaction  parameter  y , Dufour  number  Dll  Magnetic  field 
parameter  M ,Heat  absorption  parameter  (f)  , Radiation 


absorption  parameter  Q{  .Radiation  parameter  F,  Porous 

permeability  parameter  <f)x  .Solutal  Grashof  number  Gr 

mass  Grashof  num  Gift  are  presented  in  figures  1-15. 
Throughout  the  calculations  the  parametric  values  are 

chosen  as  Pr  = 0.7 1, A = 0.5, £ — 0.02, n = 0.l,t  = 1. 


Fig.l:  Velocity  Profiles  for  various  values  ofF  Fig.2:Sc  — 0.6,  Gr  = 4.0,  Gm  — 2.0,  Q1  = 2.0,  0 = 
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witli,Gr  — 4.0,  Sc  — 0.6,  Gm  — 2.0,  Qx  = 2.0,0  = 
1 ,M  = 2.0,  Du  = 0.5,0!  = 0.3, y = 0.5. 


1,  M = 2.0,  Du  = 0.3,  0!  = 0.3,  F = 2.0. 


y y 


Fig.  3:  .Velocity  profiles  for  various  values  of 
Mwith,Sc  — 0.6,  Gr  — 4.0,  Gm  — 2.0,  Q1  = 2.0,  0 = 
1.0, Du  = 0.3,0!  = 0.3,  F = 2.0,  y = 1. 


Fig.4:  Velocity  profiles  for  various  values  of 
0withSc  = 0.6,  Gr  — 4.0,  Gm  — 2.0,  Q1  = 2.0,  Du  — 
0.3,  0i  = 0.3,  F = 2.0,  y = 1,  M = 2. 


y 

Fig. 5:  Velocity  profiles  for  various  values  of 
QpwithSc  — 0.6,  ,Gr  = 4.0,  Gm  = 2.0,  M — 2.0,  Du  — 
0.1,  0i  = 0.3,  F = 2.0,  y = 1,0=  1.0. 


Fig. 6:  Velocity  profiles  for  various  values  of  Gr 
withSc  = 0.6,  Gm  = 3.0,  M = 2.0,  Qt  = 2.0,  Du  = 
O.3,0i  = 0.3,  F = 2.0, y = 1,0=  1.0 
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y 

Fig.7:  Velocity  profiles  for  various  values  of 
GmwithSc  — 0.6,,  Gr  = 4.0,  M — 2.0,  Qt  — 2.0,  Du  = 
0.3,  0!  = 0.3,  F = 2.0,y  = 1,0  = 1.0. 


Fig.  8:  Velocity  profiles  for  various  values  of 
0]  withSc  — 0.6,  Gr  — 4.0, 

Gm  = 3.0,  M = 2.0 ,Q1  = 2.0,  Du  = 0.3,  F = 2.0,  y = 
0.5,0  = 1.0. 


Fig.  9:  Velocity  profiles  for  various  values  ofSc 
withGr  = 2.0,  Gm  = 2.0,  M — 2.0,  Q1  — 2.0,  Du  = 
0.1,0!  = 0.3,  F = 2.0,y  = 0.5,0  = 1.0. 


y 


Fig. 10:  Velocity  profiles  for  various  values  Du 
withSc  = 0.6,  Q1  — 2.0,  Du  = 0.5,0!  = 0.3,  F - 
2.0,0  = 1.0,y  = 0.5,  M = 2.0. 
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Fig.ll:  Temperature  profiles  for  various  values 
ofFwithSc  = 0.6,  Gr  = 4.0,  Gm  = 3.0,  M — 2.0,  Q1  = 

2.0,  Du  = 0.5,0!  = 0.3,0  = 1.0. 


Fig.l2:Temperature  profiles  for  various  values 
of0withSc  = 0.6,  Gr  = 4.0,  Gm  = 3.0,  M — 2.0,  Q1 

2.0,  Du  = 0.5,  0!  = 0.3,  F = 2.0, 

Y = 1.0. 


Fig.  1 3 .'Temperature  profiles  for  various  values  Fig.  14:Temperature  profiles  for  various  Du  withSc  — 

ofQpwithSc  = 0.6,  Gr  = 4.0,  Gm  = 3.0,  M = 0.6,  Q1  = 2.0,  F = 2.0,  y = 0.5. 

2.0,  Du  = 0.5,  F = 2.0,0  = 1.0,  y = 0.5. 
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Fig.l  5 .'Concentration  profiles  for  various  values  of 
Sc  with  y = 0.5. 


Table. 1:  Comparison  of  present  results  with  those  of  Kim  [23]  & Pal  and  Talukdar  [6]  with  different  values  of  C ^ 


Nux/Rex 


Kim[23] 

(G  = 2.K  = 

™,UP  =0) 

Pal  and  Talukdar  [6] 

/Gr  = 2,0  = 0,  0!  = 0,F  = 0A 

V Qt  — 0,  Gm  = 0 / 

Present  Results 

Gr  = 2,  0 = 0,  0!  = 0, 

F — 0,  Q1  = 0,  Gm 
— 0 ,Du  — 0 

M Cfx  Nux/Re 

M Cfx  NuJ Rev 

M 

Cf,  Nux/  Rex 

0 4.5383 

-0.9430 

0 

4.5383  -0.9430 

0 

4.5383  -0.9430 

2 3.9234 

-0.9430 

2 

3.9234  -0.9430 

2 

3.9234  -0.9430 

5 4.4457 

-0.9430 

5 

4.4457  -0.9430 

5 

4.4457  -0.9430 

10  5.2976 

-0.9430 

10 

5.2976  -0.9430 

10 

5.2976  -0.9430 

Table:  2:  Comparison  of  present  results  with  those  of  Pal  and  Talukdar  [6]  with  different  values  of  C J,  Nllx  / Re^  (F  = 


Y Cfx  Nu  J Rcx  Shx/Rex 

C„  NuJ  Re,  ShJ  Re, 

(Pal  and  Talukdar  [6]) 

(Present  results) 

2,0  = 1,  M = 2 ,Q1  = 2 ,Gr  = 4,Gm  = 2, Sc  = 0.6,0!  = 0.3, Du  = 0.0) 
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0.0  0 

4.0441 

-1.3400 

-0.8098 

4.0442 

-1.3400 

-0.8098 

0.50 

3.7512 

-1.4825 

-1.1864 

3.7513 

-1.4825 

-1.1864 

0.75 

3.6744 

-1.5226 

-1.3178 

3.6744 

-1.5227 

-1.3178 

1.00 

3.6149 

-1.5546 

-1.4325 

3.6149 

-1.5546 

-1.4326 

Table. 3:  Effects  of  (j)f  F = 2 )and  F ( (f>\  =0.3)for  fixed  values  of  Gr  — 4,  Gm  — 2,  Sc  — 0.6,  M = 2, 0 = 1,  Q±  = 2,y  = 

0.1,  Du  = 0.5 


0i  c fx 

FCf* 

NuJ  Re, 

0.0 

6.2727 

1.0 

3.6557 

-0.8319 

0.1 

4.9698 

2.0 

3.5111 

-1.3014 

0.3 

3.5111 

3.0 

3.4129 

-1.6742 

0.5 

2.7143 

4.0 

3.3400 

-1.9860 

The  effect  of  radiation  parameterFis  shown  in  fig.  l.It  is 
observed  in  this  figure  that  velocity  profiles  decrease  with 
an  increase  in  the  radiation  parameterF.  It  is  clear  that,  the 
increase  of  the  radiation  parameter  leads  to  decrease  the 
momentum  boundary  layer  thickness  and  reduces  the  heat 
transfer  rate  in  the  presence  of  thermal  and  solutal 
buoyancy  forces.  Fig. 2 illustrates  the  effect  of  chemical 
reaction  parametery.lt  is  clear  that  at  first  the  velocity 
increases  and  later  on  decreases  uniformly  with  the 
increasing  values  ofy.  It  is  interesting  to  observe  that  the 
peak  values  of  the  velocity  profiles  attain  near  the  porous 
boundary  surface. The  velocity  profiles  for  different  values 
of  magnetic  parameter M are  depicted  in  fig.  3.  From  this 
figure  it  is  clear  that  as  the  magnetic  field  parameter 
increases,  the  Lorentz  force,  which  opposes  the  fluid  flow 
also  increases  and  leads  to  an  enhance  deceleration  of  the 
flow.  This  result  qualitatively  agrees  with  the  expectations 
since  the  magnetic  field  exerts  retarding  force  on  the  free 
convection  flow.  Fig. 4 represents  the  fluid  velocity 
decreases  as  the  heat  absorption  parameter  0 increases 
and  hence  momentum  boundary  layer  thickness  decrease. 
Figs.  5 and  13  display  the  effects  of  the  radiation 
absorption  parameter  Q}  on  velocity  and  temperature 
fields.  It  is  obvious  from  the  figures  that  an  increase  in  the 
absorption  radiation  parameter  Q1  results  an  increase  in 
the  velocity  and  temperature  profiles  within  the  boundary 
layer  as  well  as  an  increasee  in  the  momentum  and 
thermal  thickness.  This  is  because  the  large  values  of  Q1 


correspond  to  an  increased  domination  of  conduction  over 
absorption  radiation,  thereby  increasing  buoyancy  force 
and  thickness  of  the  thermal  and  momentum  boundary 
layer.  From  fig. 5,  it  is  clear  that  the  velocity  starts  from 
minimum  value  of  zero  at  the  surface  and  increases  till 
attains  the  peak  value.  For  different  values  of  the  thermal 
buoyancy  force  parameter  Grand  solutal  buoyancy  force 
parameter  Gmare  plotted  in  figs.  6 and  7.  As  seen  from 
figures  that  maximum  peak  value  is  attained  and 
minimum  peak  value  is  observed  in  the  absence  of 
buoyancy  force.  This  is  due  to  fact  that  buoyancy  force 
enhances  fluid  velocity  and  increases  the  layer  thickness 
with  increase  in  the  value  of  GrorGm. Fig. 8 illustrates  the 
variation  of  velocity  distribution  across  the  boundary  layer 
for  various  values  of  the  slip  parameter01.  It  is  observed 
that  the  velocity  increases  near  the  source  and  reaches  the 
free  stream  condition.  The  permeability  0,  is  directly 
proportional  to  square  root  of  the  actual  permeability  K. 
Hence,  an  increase  in  <f)1  will  decrease  the  resistance  of 
the  porous  medium  which  will  tend  to  accelerate  the  flow 
and  increase  the  velocity. 

The  influence  of  Schmidt  number  Sc  on  the  velocity 
profiles  are  shown  in  fig. 9.  We  observe  in  fig. 9 that  at 
very  low  values  of  Schmidt  number  (e.g.,5c  = 0.16), 
there  is  an  increase  in  the  peak  velocity  near  the 
plate  (y  — 1)  . Whereas  for  higher  values  of  Schmidt 
number,  the  peak  shifts  closer  to  the  plate.  Further  it  is 
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observed  that  the  momentum  boundary  layer  decreases 
with  an  increase  in  the  value  of.Sc. 

The  temperature  profiles  for  different  values  of  the 
radiation  parameter  F are  shown  in  fig.ll.lt  is  observed 
that  an  increase  in  the  radiation  parameter  F results  a 
decrease  of  the  thermal  boundary  layer  thickness.  Further 
it  is  observed  from  that,  the  temperature  is  very  high  at  the 
porous  boundary  and  asymptotically  decreases  to  zero 
asy->oo.  Fig.  12  depicts  the  variations  in  temperature 
profile  against  spanwise  co-ordinate  y for  different  values 
of  heat  absorption  parameter  0 . From  this  figure,  it  is 
clearly  understood  that  heat  absorption  parameter 
condenses  the  thickness  of  the  temperature  boundary  layer 
because  when  heat  is  absorbed,  the  buoyancy  force 
decreases  the  temperature.  The  effect  of  Schmidt  number 
Sc  on  the  species  concentration  profiles  is  shown  in 
fig. 15. It  is  clear  that  the  concentration  decreases 
exponentially  and  attains  free  stream  condition.  Also  it  is 
noticed  that  the  concentration  boundary  layer  thickness 
decreases  with. Sc. 

The  temperature  profiles  for  different  values  of  Dufour 
number  in  heat  absorption  (0  > 0)  and  heat  generation 
(0  < 0)  cases  are  shown  in  fig.  14.  It  can  be  seen  that  the 
fluid  temperature  increases  with  Dufour  number  in  both 
the  cases.  Thermal  boundary  layer  thickness  is  higher  in 
the  case  of  heat  absorption  than  that  in  the  case  of  heat 
generation.  Physically,  the  Dufour  term  that  appears  in  the 
temperature  equation  measures  the  contribution  of 
concentration  gradient  to  thermal  energy  flux  in  the  flow 
domain.  It  has  a vital  role  in  enhancing  the  flow  velocity 
and  the  ability  to  increase  the  thermal  energy  in  the 
boundary  layer.  As  a result,  the  temperature  profile  at  all 
time  stages  increases  with  the  increase  in  Du. 

Fig.  10  displays  the  velocity  profiles  for  various  values 
Dufour  number  in  the  cases  of  cooling  and  heating  of  the 
plate.  The  velocity  increases  with  Dufour  number  in  the 
case  of  heating  of  the  plate  and  opposite  trend  is  observed 
in  the  case  of  cooling  of  the  plate. 

From  Table.  1 and  Table. 2,  it  is  observed  that  the  skin 
friction  at  the  plate  decreases  with  increasing  chemical 
reaction  parameter  y or  Magnetic  parameterM.  It  is  also 
noticed  that  the  Nusselt  number  and  Sherwood  number 
decreases  with  increase  in  y. 

The  influence  of  slip  parameter  0X  on  skin-friction,  and 
radiation  parameterFon  skin-friction  and  Nusselt  number 

Appendix 

p _ Sc+J(_Sc2+4Scy'j  _ Pr+/Pr2  +4(F+0)Pr  _ l+74(M1+n) 


are  presented  in  Table.  3.  It  is  observed  that  the  skin 
friction  at  the  plate  decreases  with  increasing  0,  or  F.  As 
the  radiation  parameter  increases  the  heat  is  transfered 
from  the  plate  to  the  fluid. 

V.  COMPARISON  OF  RESULTS 

In  order  to  examine  the  accuracy  of  the  results  of  the 
present  study,  it  is  considered  that  the  analytical  solutions 
obtained  by  Kim  [23],  Pal  and  Talukdar[6]  who  computed 
the  numerical  results 

for  skin  friction  coefficient  and  Nusselt  number.  These 
computed  and  compared  results  are  presented  in  Table.  1. 
From  this  table  it  is  interesting  to  observe  that  the  present 
results  in  the  absence  of  Diffusion-thermo  effect  and 
Porous  medium  are  in  good  agreement  with  the 
corresponding  results  obtained  from  Pal  and  Talukdar  [6]. 
It  is  also  observed  that  the  present  results  are  in  good 
agreement  with  those  of  Kim  [23]  when  G — 2,  up  — 
0,  Du  — 0,  Gm  — 0.  Also  from  Table.  2,  it  is  clear  that  the 
present  results  for  skin-friction,  Nusselt  number  and 
Sherwood  number  for  different  values  of  chemical 
reaction  in  the  absence  of  porous  medium  and  diffusion 
thermo  effect  are  in  good  agreement  with  the 
corresponding  results  of  Pal  and  Talukdar  [6],  which 
clearly  shows  the  correctness  of  our  present  analytical 
solutions  and  computed  results. 

VI.  CONCLUSIONS 

In  this  paper  we  have  studied  the  Dufour  effect  on  an 
unsteady  MHD  convective  heat  and  mass  transfer  flow 
through  a high  porous  medium  over  a vertical  porous 
plate.  From  the  present  study  the  following  conclusions 
can  be  drawn. 

1.  The  diffusion-thermo  parameter  increases  the 
thermal  and  momentum  boundary  layer  thickness. 

2.  The  fluid  velocity  increases  with  an  increasing 
values  of  slip  parameter01. 

3.  The  skin  friction  at  the  plate  increases  with  an 
increasing  values  of  Du. 

4.  The  rate  of  heat  transfer  coefficient  at  the  plate 
(Nux)  increase  with  Dufour  number  while  it  is 

decreases  as  chemical  reaction  parameter  (y)  or 
Radiation  parameter  (F)  increases. 
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m, 


_Sc+]Sc2  +4 Sc(y+n)  ___  _ Pr+1/f’r2+4Pr(F+0+n)  _ _ 1+7 1+4 M1  A _ (- PrQ1-DuPl 2) 

L — ,?72r  — •Ttifi  — , /li  — o " " , 

b 2 b 2 0 2 1 (P^-PrP!-(F+0)Pr) 


. (Gr^x+Gm)  , _ G7~(l-i4x)  >1  _ 

2 _ (P^-Pi-MO’^3  _ (m2-m2-M1)  ’ 4_ 


Am^B^ 


A*=- 


AScP1 


(Pz  -Px.Sc-Sc(y+77)) 


A _ (-PrthAs-DuAsPj+PrPiAAi)  _ 
6 (P^— PrPi— Pr(F+0+n))  ’ 7 


(. PrAm2-PrAA1m2 ) 
(m2-m2Pr-Pr(F+0+n)) 


'>  ^8  — 


t-PrQ1+PrA5Q1-Dum%+DuA5mA') 
(mJ-m4P? — Pr(F+0+n)) 


i49  — 1 (^46  + i47  + i4g),i410  — 
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7774  -7774  - (Mx  +7i) 


(1  + 0x7776  ) 


B,  = 


(l+01m3 ) 


£1  = 
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